home *** CD-ROM | disk | FTP | other *** search
-
-
-
- ____BBBBAAAAKKKKVVVVEEEECCCC((((3333FFFF)))) ____BBBBAAAAKKKKVVVVEEEECCCC((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- BAKVEC, SBAKVEC - EISPACK routine. This subroutine forms the
- eigenvectors of a NONSYMMETRIC TRIDIAGONAL matrix by back transforming
- those of the corresponding symmetric matrix determined by FIGI.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee bbbbaaaakkkkvvvveeeecccc((((nnnnmmmm,,,, nnnn,,,, tttt,,,, eeee,,,, mmmm,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr mmmm,,,, nnnn,,,, nnnnmmmm,,,, iiiieeeerrrrrrrr
- ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn tttt((((nnnnmmmm,,,,3333)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,mmmm))))
-
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssbbbbaaaakkkkvvvveeeecccc((((nnnnmmmm,,,, nnnn,,,, tttt,,,, eeee,,,, mmmm,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr mmmm,,,, nnnn,,,, nnnnmmmm,,,, iiiieeeerrrrrrrr
- rrrreeeeaaaallll tttt((((nnnnmmmm,,,,3333)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,, mmmm))))
-
-
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- On INPUT
-
- NNNNMMMM must be set to the row dimension of two-dimensional array parameters
- as declared in the calling program dimension statement.
-
- NNNN is the order of the matrix.
-
- TTTT contains the nonsymmetric matrix. Its subdiagonal is stored in the
- last N-1 positions of the first column, its diagonal in the N positions
- of the second column, and its superdiagonal in the first N-1 positions of
- the third column. T(1,1) and T(N,3) are arbitrary.
-
- EEEE contains the subdiagonal elements of the symmetric matrix in its last
- N-1 positions. E(1) is arbitrary.
-
- MMMM is the number of eigenvectors to be back transformed.
-
- ZZZZ contains the eigenvectors to be back transformed in its first M
- columns. On OUTPUT
-
- TTTT is unaltered.
-
- EEEE is Destroyed.
-
- ZZZZ contains the transformed eigenvectors in its first M columns.
-
- IIIIEEEERRRRRRRR is set to Zero for normal return, 2*N+I if E(I) is zero
- with T(I,1) or T(I-1,3) non-zero.
- In this case, the symmetric matrix is not similar
- to the original matrix, and the eigenvectors
- cannot be found by this program. Questions and comments should be
- directed to B. S. Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL
- LABORATORY
-
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-